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Solutions of nonlinear multi-component Euler-Monge partial differential equations are constructed in n spatial 
dimensions by dimension-doubling, a method that completely linearizes the problem. Nonlocal structures are an 
essential feature of the method. The Euler-Monge equations may be interpreted as a boundary theory arising 
from a linearized bulk system such that all boundary solutions follow from simple limits of those for the bulk. 



For any theory^ with an infinite number of conservation laws, we may always assemble the conserved currents 
into a generating function involving a spectral parameter a. If that spectral parameter is independent of any other 
spacetime dimensions in the theory, as is possible in the simplest cases, then effectively the theory possesses an 
extra dimension^. Moreover, it is always possible to openly include this extra dimension in some of the dynamical 
OO ' equations, and not just leave it as a dimension sub rosa. 

For example, suppose a theory is originally expressed in terms of coordinates (x, t) with an infinite number of 
conserved currents: d t p^ (x,t) — d x J^ {x,t), n G N. Then by defining p(x, t, a) = J2 n ( n + l) a "/^™' 1 (x,t), as 
opposed to J2n a 'V'"' 1 { x : t), and J [x, t, a) = J2 n a ' l+1 ^™' 1 (x, i), as opposed to J2n o- n J^ i x i t), we have rendered all 
' the conservation laws as a single second-order higher-dimensional partial differential equation (PDE): dtp(x,t,a) = 
d x d a J [x, t, a), as opposed to the first-order dtp (x, t, a) — d x J (x, t, a). Hence our choice for the current generating 
functions has fully exposed an extra dimension in the PDEs satisfied by those generating functions. The extra 
dimension here does not just ride along as a suppressible label for the currents but it appears explicitly, perhaps 
even unavoidably, in the dynamical equations. Of course this immediately raises issues about whether the theory 
requires a to appear explicitly for all dynamical equations to be cogently expressed in terms of the original plus extra 
dimensions, and about covariance properties for the theory in the complete set of dimensions. 

In this paper we address these issues for a simple but very generally applicable class of nonlinear PDE's [l7[ |: 
The first order Euler-Monge (E-M) equations <9 t u = (u • V) u. We find the full dynamics of these nonlinear theories 
are elegantly encoded into a higher dimensional set of linear "heat" equations obtained through dimension doubling 
(x) — > (x, a), where for each spatial coordinate Xi there is an associated coordinate given by spectral parameter a^. 
The original dynamical variables are obtained as spectral parameter boundary limits, Ui (x, t, a) a _ >0 ' m (x, t). The 
fact that the higher dimensional theory is linearized strongly argues that this is the right approach to take. In the 
linearized theory, the pairs (xi, ai) act like "light-cone" variables in the enlarged set of dimensions such that the heat 

equations for all the dynamical variables are of the form (d/dt — 2^=i d 2 / dojdx^j Ui (x, t, a) = 0. Thus the extra 
dimensions appear explicitly and, indeed, unavoidably in these linearized dynamical equations. 

We also find Nambu brackets 18 of the fields, of all orders up to the full Jacobian, as a remarkable feature of the 
linearizing maps. We know of only one other field theoretic example || where these brackets appear so naturally. 
Moreover, the linearizing maps are nonlocal in all but the simplest, one component case. The nonlocal structures 
appropriate for E-M equations with two components in two spatial dimensions are evocative of phase factors in Wilson 
loops (cf. strings), and when the E-M equations describe n component fields in the original n spatial dimensions 
these structures extend to higher dimensional constructions involving integrals over n — 1 dimensional submanifolds 
(cf. (n — l)-branes). In the one dimensional, one component case, the E-M solution is obtained algebraically from 
the dimensionally-doubled "bulk" solution for all values of the single spectral parameter. In higher dimensional or 
multi-component cases the dependence of the solutions on the spectral parameters is more involved. Nevertheless, 
in all cases the solutions of the E-M equations may be obtained from simple limits of those for the bulk. 

1 We dedicate this paper to Peter Freund on the occasion of his becoming Professor Emeritus at the University of Chicago, and thereby 
begin with general remarks about the origin of extra dimensions, allowing for the possibility that these are similar to but not necessaril 
on the exact same footing as the original dimensions. For related points of view within the Kaluza-Klein physics framework Jlj], see 

2 More precisely, an extra bosonic dimension. A finite number of conservation laws evokes extra fermionic or anyonic dimensions. 
involving fcth order superspace or Grassman variables. This and other non-commutative geometries will not be discussed further here 
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The Euler-Monge equations first appeared in 18th and 19th century studies of fluid dynamics |T(J and analytic 
geometry jlT]] . Riemann took up a study of the equations in the context of gas dynamics, discussing the equations as 
a theory of invariants pp[ (for a modern textbook treatment, see ||). His approach is widely applicable to almost all 
nonlinear flow problems, although it does not triumph over turbulence. A systematic modern discussion of the E-M 
equations that synthesizes ideas from both geometry and invariance theory can be found in the review by Dubrovin 
and Novikov . Most contemporary texts and reviews stress the universal role played by these nonlinear transport 
equations in accordance with Whitham's theory [ pT| . Essentially all nonlinear waves, even those in dispersive and 
dissipative media, involve E-M equations, or simple variants of them, if the nonlinear wavctrains are slowly varying. 
This makes the equations particularly useful for analyzing the asymptotic behavior of nonlinear solutions. The E-M 
equations and their conservation laws also serve as a useful starting point in Polyakov's study of turbulence jl^] but 
without yet leading to a general solution of the Navier-Stokes equations. 

The first order E-M equation ^ = it|j also gives rise to the Bateman equation upon substituting u — §^j§^- 

The resulting second order nonlinear PDE is 4> 2 x (t>u ~ < ^ ( t ) x ( t > t ( i > tx + ^t^xx = ^ ano - ^ s weu known to possess a general 
implicit solution given by solving tSo((f>) + xSi(4>) — constant, where So and Si are arbitrary differentiable functions 
of 4>(x,t). The structure of this solution incorporates the covariance properties of the PDE: If ^ is a solution, so 
is any function of <p. In fact, curiously, the generalization of this solution to n + 1 functions Sq(^), Si(<p) of <fi (x, t), 
x = [x%, • ■ • ,x n ), subject to a single constraint tSo(tfi) +2 x i ^i( t t > ) = 0, is a "universal solution" |14| to any equation 
derived from a Lagrangian which is homogeneous of weight one in the first derivatives of (f>. 

Thus the Euler-Monge equations appear widespread across a very broad landscape of physics and applied math- 
ematics problems, and therefore it is important to understand their solutions at as many levels as possible. To that 
end we shall map all solutions of the E-M equations in arbitrary dimensions into solutions of second-order linear 
equations. This type of map is reminiscent of the Cole-Hopf J(], ^| transformation (thoroughly reviewed in fl6|| ) used 
to linearize the Burgers ^, ^| nonlinear diffusion equation, but there are important differences here. The Cole-Hopf 
transformation only works for curl- free u, does not use extra dimensions, and fails for = n (the diffusivity) . The 
map to follow works for all u, curl-free or otherwise, does use extra dimensions, but works only for k = 0. (We hope 
to extend the method to k ^ in a subsequent study.) 

We believe it is most efficient to present our results summarily as a small set of Theorems, for which we sketch 
the salient features of their proofs. In the following, Ai n is the n dimensional nonlinear Euler-Monge operator and 
TL n is an associated hyperbolic heat operator (introduced in Jl9|]). 

9 d „, d A d 2 



dt dxj dt ^— ' dxjdaj 

j— 1 3 — 1 

To begin, however, we will generalize these two definitions to allow for an arbitrary function F in the most elementary 
results in one spatial dimension. We fincf] 

Theorem 1 in 1: (x,t, a) — F (^) J^U {x,t, a) if and only if j^u (x,t) = F (u (x,t)) J^u (x,t) where 

e au(x,t) _ ^ ^ 

U {x, i, a) = , u{x,t) — — In (1 + aU (x, t, a)) 

a a 

and F is any function with a formal power series. 
Proof of 1 in 1: By direct calculation 

(I AD s) ^ — (!»<*■«>- *<■"-<» £•<*<>) 

and the Theorem follows. ■ 

Corollary 1 in 1: The formal solution for U (x, t, a) in terms of U (x,t — 0, a) is given by 

f e -n(x,t) _ i) / a = e *^(£)£ (Uu{x) _ A /a \ 

with u (x, t = 0) = u (x) . This is an elementary consequence of the Theorem. ■ 

3 We enumerate the Theorems by "k in d" , where k is the number of field components and d is the number of spatial dimensions in 
the Euler-Monge equations. 
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The bulk solution U (x, t, a) may also be viewed as a simple one-parameter deformation of the boundary data 
u (x, t), with the extra dimension serving as the deformation parameter. In this exceptional one-component case, we 
may easily extract u (x, t) from U (x, t, a) for any value of the extra dimension a as given by the logarithmic expression 
in the Theorem. But in particular, we may extract u(x,t) as a limit of the bulk solution u(x,t) = lim U (x,t,a). 

This immediately yields the time series solution |u| to the previous E-M equation as the limit 



, +N t F(JL)JL fe au ^ -l\ 

, (x, t) = Um e v »» / a* 

q^o \ a J 



3=0 v Jl 



where we assume F (locally) invertible in the last stepf^]. Similar Corollaries and time series solutions are obvious 
consequences of all our results, and may be incorporated directly into each Theorem. For example, one independent 
field u in spatial dimensions [x, y\, ■ ■ • , y n ) with dependent "velocity fields" (u, v\ (u) , ■ • ■ ,v n (u)) leads to 

Theorem 1 in (n+1): (x, y, t) = u (x, y, t) -§^u (x, y, t) + v i ( u ( x > Y, *)) ( x > Yj t) if and onl y if 

/■«(*.y.*) / " \ t( * +t- 32 ) r u(x ' y) ( A \ 

J d ue XP Uu + ^ bM u) =e^ ^* =1 ^1^1) J du exp I au + > bjVj (u) J 

Proof of 1 in (n+1): By direct calculation, with U (x, y, t, a, b) = J^ x ' y '^ du e xp (au + Y^7=i ^ iVi ( u )) > 
( d d 2 " d 2 \ 

U"^ _ ^%^r (:E ' y ' < ' a ' b)= 



Id d 9 \ I \ 

[~di U ^' Y ' ^ ~ " ^' y ' ^ fe U ^' y ' ^ ~ ^ ^ ^ ^' y ' ^ % W ^' y ' ^ J 6XP ( ° U ^' y ' ^ + ^ ^ ^ ^' y ' *^ J 

So, as given, the higher dimensional heat equation is satisfied by the integral form U (x, y, t, a, b) if and only if the 
given one-component generalization of the E-M equations holds. The RHS of the relation in the Theorem is then 
just the formal solution of the heat equation, as in the previous Corollary]^ ■ 

The last result does not allow for a simple extraction of u(x,y,t) from the integral form of U (x, y, t, a, b) for 
non- vanishing a, b. However, it does have the simple limit lim U (x, y, t, a, b) = u (x, y, t), so extraction is trivial 

on the boundary a, b — > 0. This is true of all the heat equation solutions to follow. Also note, U (x, y, t, a, b) in this 
one-component case is an integral over the field value. Nevertheless U is still local in all the dimensions, no matter 
how many. Locality in the original spatial dimensions will not hold, however, for maps of multi-component fields in 
higher dimensions. This is first illustrated by 

Theorem 2 in 2: H 2 U = H 2 V = if and only if M 2 u = M 2 v = where (e(s) = ±\ for s ^ ) 

U (x, y, t, a, b) = dre(y- r) e «»(*.r,t)+K*,r,t) du{x,r,t) 

V (x, y, t, a, b) = J™ dqe(x- q) e ^( q ,y,t)+bv{ q ,y,t) MMl_ 



4 Given the close relation of the Monge and Bateman equations, it might be expected that the latter also admits a power series solution 
of simple form. Indeed this is so. Treating the equation as hyperbolic with initial conditions </)(a;,0) = f(x), d(j>(x,0)/dt = g(x), the 
time series solution of the Bateman equation is 



_■ i t J- ' 



£i [l + 3)ldxi \ {df(x)/dxy 

16] 

obtained by elimination of u from the first order equations jlj] 



As is true for the Bateman equation and the one-component Monge equation in one spatial dimension, there is a corresponding 
second order equation for the 1 in (n-)-l) case which our soiutipji satisfies. It is the so-called "Universal Field Equation" which may be 
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Proof of 2 in 2: Again by direct calculation, assuming u, v, and their derivatives vanish asymptotically in x, y, 

H 2 U(x,y,t,a,b) = e au ( x ^ +bv(x ^M 2 u(x,y,t) 

+b J°° dre(y- r) e ««(*.r.*)+K*.r,t) ( ^ U ^ l) M 2 v (x, r, t) - — M 2 u (x, r, t) 

n 2 V(x,y,t,a,b) = e ^(x, v ,t)+bv(x, y ,t) M2V ^ ^ t j 

+a f dqe{x- q) e»<™.0+*-(™,*) ( ^^W (q, y, t) - M 2 v (q, y, t) 

The complete Theorem then follows by also using the obvious limits lim H 2 U (x,y,t,a,b) — M 2 u(x,y,t) and 

a.b^O 

lim 7i 2 V(x,y, t, a, b) = M 2 v(x,y, t). ■ 

As advertised, the two-component map in two spatial dimensions involves a nonlocal transformation between E-M 
and heat equation solutions: It features line integrals over the original spatial variables. The map is still local in 
the extra dimensions, however. This nonlocality in the original dimensions persists and is even extended when more 
components and more spatial dimensions are considered. As a further illustration before giving the generalization 
to an arbitrary number of dimensions, we have 

Theorem 3 in 3: H 3 U = H 3 V = H 3 W = if and only if M 3 u = M 3 v = M 3 w = where 

U (x, y, z, t, a, b,c) = J dr e (y - r) e au+bv+cw du( - x £ z >^ _ C JJ drds £ ( y - r ) e (z - s) e au+hv+cw {«, w} rs (x, r, a, t) 
V (x, y, z, t, a, b, c) = Jdse(z- s) e au+bv+cw dv ^ s ^ - a JJ dqds £ {x-q)e{z-s) e au+bv+cw {v, u} sq (q, y, s, t) 
W(x,y,z,t,a,b,c) = jdqe(x- q) e au+bv+cw — ~ b J J dqdr e (x - q) e (y - r) e au+bv+cw {w, v} qr (q, r,z,t) 

Proof of 3 in 3: There are a few essential new ingredients needed to complete the proof by direct calculation 
in this case. Define Poisson brackets as usual by 

, , dudv du dv 

i U > V ir S = g^-g; ~ 

where u and v are any two functions of the independent variables r and s. Then it is straightforward to show 

§t {w ' v} *y ~ Tx ( u {m ' v}z y) "k( v {u ' v}z y) ~§~z( w {m ' v} *y) = {M3U ' v} *y + {u ' Msv} *y 

Ft {u ' v] *y ~lL( u {u ' v} *y) ^Fj( v {u ' v} *y) ~§~z( w {m ' v} *y) = {Mu > + {u ' Mv} *y ~ {w ' v ' w] ^ z 



as well as similar relations obtained by permutation of dependent and independent variables. In the last relation 
we have introduced the totally antisymmetric Nambu triple bracket (i.e. Jacobian, in this 3-dimensional case) 

du du du du dv dvj 

{«> v, w} xyz = — {v, w} yz + — {v, w} zx + — {v, w} xy = — {v, w} yz + — {w, u} yz + — {«, v} yz 

Once equipped with such relations, the complete proof of the Theorem is tedious, perhaps, but not subtle. (See the 
generalization to follow for additional details.) ■ 

The nonlocality appearing in our map for three components in three spatial dimensions is two-dimensional: It 
features surface integrals over pairs of the original spatial dimensions, perhaps evocative of membrane-based phase 
factors. Nonetheless, the map is still local in the extra dimensions and the E-M solutions are again trivially given 
by boundary limits of the bulk constructions. The nonlocality is extended to (n — l)-dimcnsional integrals when 
n-component linearizing maps are constructed in n spatial dimensions. This is explicit in 
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Theorem n in n: H n U k (x, t, a) = if and only if M n Ui (x, t) = for i, fc e {1, • • • , n} where 



U k (x, t, a) 



dgi • ■■dq n 5(q k - x fc ) 



x det 



V ( e x„)e a ""' 1 ) ••• ^ (e(q n x n )e a " Un ) j cxc i udo feth row 

and /cth column 



Proof of n in n:[] Consider only the first component (et sic de similibus). 



tfi(x,i,a) 



dqi ■■ - dq n 5 (qi - Xi) 



«i 



det 



£-(e (q n -x n )e a ^) 



^(e (q 2 -x 2 )e a ^) 



^(e (q n -x n )e a ^) 



dq 2 --- dq n 



dq 2 --- dq n 



gOlUl _ i 

ai 

gaiiii _ ^ 

ai 



()</ 



S (e(q 2 -x 2 )e a2U *) 



9 (e(q 2 -x 2 )e a ^) 



det 



(e (q n -x n )e a ^) 



(e ((7„-x„)e a " M ") 



dq. 



12 



dq, 



(xi,q 2 , ■■■ ,q n ,t) 



(e (<7„-x„)e a " u ") 



»3 



where in the last expression the i k dummy indices, k S {2, • • ■ , n}, are summed from 2 to n, i.e. 1 is excluded. Now 
we integrate by parts assuming all fields and their derivatives vanish as x — ► 00. To do this, there are clearly n — 1 
equivalent choices. We elect to integrate t^- by parts to obtain 



Ui (x, t, a) 



9 

d<72 • • • dqn £ (<& ~ x 2 ) e a2 " 2 - — 



dq 2 ■ ■ ■ dq n x 



1\ d 



,^-(s(q 3 -x 3 )e a ^)--- — 
ai J dq l3 dqi 



xe (q 2 - x 2 ) ^— - ( <5 i3 3^ (33 - 2:3) + a3£ (93 ~ x 3 ) 

dq l2 \ dq n 



Si nn S (q n ~x n ) + a n e (q n - x n ) 



du 



dq 



n \ e a-u 



Expanding out the products of the various paired terms in parentheses in the last line gives 
Ui (x, i, a) = —Ei 2 ...i n a 3 ■ ■ ■ a n 



, , . / \ / \ / duxdu-i du„ au 
dq 2 dq 3 ■ ■ ■ dq n e [q 2 - x 2 ) e [q 3 - x 3 ) ■ ■ -e (q n - x n ) - — — e 



dq i2 dq i3 dqi 



( 



'(I p 71 / p 

£ i2---i n dq 2 e (q 2 - X 2 ) Sjij Y\ ( a fc / d 1k £ (fl 

j=3 J \ fc=3 ^ 



k - x k ) 



du k \ 1 dux au 



\kjtj 



dqi k J dqi 



£ i2---i n ^2i^2i I dq 2 s (q 2 — x 2 ) Sjij 8ki k Yl ( 0,7,1 / dq m £ {q 

J =3fe=4- y I m =3 ^ ° qi ™/ / 

k>j \mjtj,k / 

0j y ^ 2 e (« - x 2 ) y ^ e - x 3 ) 



du m \ j dui ga . u 



9%: 



| dq 2 e(q 2 -x 2 ) |^e a u 



6 After obtaining this result, we learned that Polyakov had found inhomogeneous versions of the same higher-dimensional "heat" 
equations are obeyed by certain correlation functions in his statistical treatment of turbulence |19| (especially p 6188). 



5 



That is to say, the result is given in terms of Nambu brackets [[l8| of all ranks from n—1 down to 2 (i.e. Poisson), 
as well as a final single derivative term. Thus 

Ui (x, t, a) = -a 3 ■ • • a n J "J dq 2 dq 3 ■ ■ ■ dq n e (q 2 - x 2 ) e {q 3 - x 3 ) ■ ■ ■ e (q n - x n ) {u\,u 3 , ■ ■ • ,U n } 23 .,. n e au 
~lLi J d l2e(q 2 -x 2 ) f[ ( a k J dq k e{q k -x k ) S j {ui, u 3 , ■ ■ ■ Uj-i, Uj+i, ■ ■ ■ , u n } 23 ... j _ lj+1 ... n e a ' u 

^2^2 a i ak / dq2 e ^ 2 _ 353 ) / ^ e ^ ~ ^ / ^ e _ Xfe ) {"i' u fe}2jfe ea u 

j=3 k=i J J J 

k>j 



^2 a i J d l2 £ (q 2 - x 2 ) J dqj e (q-j - Xj) {u\,Uj} 2j e au - J dq 2 e (q 2 - x 2 ) e au 
j— 3 



In the preceding equation, it is to be understood that the sum in the second RHS row begins at its lower limit with 
-<z 4 • • • a n J ■■ ■ J dq 2 dq4 ■ ■ ■ dq n e (q 2 - x 2 ) e (? 4 - 2:4) • • • e (q n - x n ) {m, u 4 , • • • , ii„} 24 ...„ e a ' u 
and terminates at its upper limit with 

-a 3 ■ ■ ■ a n _i J ■■■ J dq 2 dq 3 ■ ■ -dq n -i e (q 2 - x 2 ) e (q 3 - x 3 ) ■ ■ ■ e (q n -i - x n -i) {ui, u 3 , ■ ■ ■ ,u n -i} 23 _.. n _ 1 e au . 

Next we act with the heat operator on Ui (x, t, a). The e's permit the appropriate "outside" (i.e. x) partials to 
be converted, through integration by parts, into "inside" (i.e. q) partials. Also, factors of ai outside the exponentials 
produce some extra terms from the cross-partials g^rg~ m Hn- We obtain 



K n Ui (x, t, a) = -a 3 ■ ■ ■ a n J ■■■ J dq 2 dq 3 ■■■dq n e (q 2 - x 2 ) e (q 3 - x 3 ) ■ ■ ■ e (q n - x n ) H n ({ui, u 3 , ■ ■ ■ , M„} 23 ... n e au ) 
^ 7^- (a 3 • • - a„) / • • • / dq 2 dq 3 ■ ■ ■ dq n — [e (q 2 - x 2 ) e (q 3 - x 3 ) ■ ■ ■ £ (q n - x n )] {u x ,u 3 , ■ ■ ■ ,Un} 23 ...„ e a u 
fn.. A 



i=3 



j=3 



\ k=3 



71 f n d d n { f 

J2 / dc t2 £ - X 2 ) g-gr II ( ak / dqk £ ~ X *) 

3=3 J i=3 1 1 \ k=3 ^ J 



{ui,V, 3 ,---Uj-i,V,j 



j+ir ■ ■ , u n} 23 ...j_ij + i... r 



Y J d< l2 £ {q 2 - x 2 ) JJ ^a k J dq k e(q k -x k ) S j^j U n ({u u u 3 , ■ ■ ■ Uj_ u u j+ i, ■ ■ ■ , U n } 23 „ d _ lj+1 

( n . , . . \ 

/ 

^ aj J dq 2 e (q 2 - x 2 ) J dqj £ (q-j - Xj) TL n {{uuUj} 2 . e au ^) 

^^"^ fe aj / dq2£{ - q2 ~ x ^ j dq j £ ( q 3~ x 3)j { u i> u ihj ea u 



i=3 
d d 



i=3 



J dq 2 e(q 2 -x 2 ) H n (|^e a u 



The first RHS line of H n Ui reduces to terms linear in the E-M equations for the m's. The second and third RHS 
lines combine to give similar terms linear in the E-M equations. And so it goes with subsequent pairs of RHS lines, 
until finally the last two RHS lines combine to give terms linear in the E-M equations. 

To establish these statements, one needs to use several identities involving the action of the heat operator on 
exponentially weighted derivatives of the component fields, in particular on so- weighted Nambu brackets. For 
example, these identities range from the simplest for the full Jacobian 



U n (e a u {wi,u 2 , • ■ • ,u„} 12 ...J = e a u (a-A4„u) {u 1} u 2l ■ ■ ■ , w„} 12 ...„ 

+ e au ({M n ui,u 2 , ■■■ , u n } 12 ... n + {ui, M n u 2 , ■■■ , u n } 12 ... n H h {u x , u 2) ■ ■ • ,M 
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to those involving lower rank Nambu brackets such as 

Tin (e aU {u 2 , « 3 , ■ • • , M„} 23 ...J = 6 aU (- {ui,U 2 , U 3 ,--- , U„} 123 ...„ + (a-X„u) {u 2 , tt 3 , ■ • • , U n } 23 ... n ) 

n w n}23-n) 

including that needed to deal with the first RHS line of H n Ui 

Tin (e a u {ui,ti3, • • • , Un} 23 ...„) = e a u (a-M n u) {ui, u 3 , ■ ■ ■ , u n } 23 ...„ 

+ e au ({M n ui,u 3 , ■■■ , Un} 23 ... n + {ui : M n u 3l ■ • ■ , u n } 23 ... n H h {ui, u 3 , • • ■ , A4 n w n} 23 ...„) 

as well as other relations obtained by permutations of the indices of these, etc., all the way down to the final 

H n {e a u d jUk (x, t)) = e a u ^dj (M n u k ) + (a-M n u) d 3 u k - W>M.,iJ 

as needed to deal with the last two RHS lines in Ji n U\. All such identities are straightforward to substantiate by 
direct calculation. 

Thus, given the E-M equations for the u's, the heat equation for U\ follows. Moreover, the only terms on the 
RHS of Ji n XJ\ which survive in the limit of vanishing spectral parameters are the last two lines, which give 

lira H n Ui (x, t, a) = M. n u\ (x, t) 

Thus, given the heat equation for U\, the E-M equation for u\ follows. Similar results obtain for all the other 
components, so that H n Uk = iff M n Uj =0. ■ 

Corollary I of n in n: Formally, time evolution in the bulk is given by 

U(x,i,a) = e t5 ^ =1 ln ^JV (^,t = 0,a) 
This gives a time-series solution on the boundary upon taking the limit a > 0. 

u(x,t) = lira e'^^ 1 ^U(x,t = 0,a) 
with initial boundary data u (x) = lim U (x, t = 0, a). 

a^O 

Corollary II of n in n: The n-fold infinite sequences of conservation laws for the E-M equations in n spatial 
dimensions are directly encoded into the bulk solutions. 

d 

— U k (x,i,a) = V ■ Jfe (x,i,a) , Jfe (x, t, a) = V a C4 (x, t, a) , k G {1,2, ••• ,n} . 
at 

Explicit sequences of charge and current densities on the boundary follow immediately from power series expansions 
in the dj. Both Corollaries are elementary consequences of the Theorem and the heat equation obeyed by U. ■ 

This is as far as we have completed the application of the higher dimensional approach to classical nonlinear PDEs. 
It remains to apply this approach to other types of nonlinear PDEs, in particular to those higher-order extensions 
of the E-M equations involving dispersion, such as the Korteweg-deVries equation, and to those involving diffusion, 
such as the Burgers and Navier-Stokes equations. Another immediately obvious challenge is to carry the method 
over to quantum field theories (QFTs). This will not be done here. However, we suspect that the implementation 
of these ideas in QFT will involve the use of quantum Nambu brackets (QNBs), given that the classical versions 
of these appear above. QNBs have a long-standing notoriety, but recently it has been shown that theirs is an 
undeserved bad reputation. QNBs can be defined in terms of operators (or in terms of non-commutative geometry) 
so as to fulfil their expected roles in the quantum evolution of dynamical systems. Perhaps these developments will 
be useful to meet the challenge of quantizing the E-M equations as well as their higher-order generalizations. 

As emphasized previously, the Euler-Monge equations appear widespread throughout physics and the mathematics 
of nonlinear partial differential equations. Based on the maps we have presented to linearize these equations, we have 
come to the following remarkable conclusion. Extra dimensions and nonlocal structures are probably universally 
applicable features to be found upon analyzing solutions of nonlinear partial differential equations, and hence quite 
natural constructs in almost all physical theories. 
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